Gravito- magnetic amplification in cosmology 



Christos G. Tsagas 
Section of Astrophysics, Astronomy and Mechanics, Department of Physics 
Aristotle University of Thessaloniki, Thessaloniki 54124, Greece 



Abstract 

Magnetic fields interact with gravitational waves in various ways. We consider the cou- 
pling between the Weyl and the Maxwell fields in cosmology and study the effects of the 
former on the latter. The approach is fully analytical and the results are gauge-invariant. 
We show that the nature and the outcome of the gravito-magnetic interaction depends on 
the electric properties of the cosmic medium. When the conductivity is high, gravitational 
waves reduce the standard (adiabatic) decay rate of the B-field, leading to its superadiabatic 
amplification. In poorly conductive environments, on the other hand, Weyl-curvature distor- 
tions can result into the resonant amplification of large-scale cosmological magnetic fields. 
Driven by the gravitational waves, these i?-fields oscillate with an amplitude that is found to 
diverge when the wavelengths of the two sources coincide. We present technical and physical 
aspects of the gravito-magnetic interaction and discuss its potential implications. 

PACS numbers: 98.80.Jk, 04.30.-w, 98.62.En 

1 Introduction 

The origin of the large-scale magnetic fields seen in the universe today remains an enigma, despite 
the efforts of the scientific community to resolve it [Tl[2]. The structure of the galactic i?-fields, 
primarily of those found in spiral galaxies, seems to support the dynamo idea fS]. However, the 
galactic dynamo requires an initial seed-field to operate and the origin of these magnetic seeds 
is still elusive. Provided the amplification is efficient, the strength of the seed can be as low as 
~ 10~^^ G by the time the host galaxy is formed [4]. This requirement could be relaxed down 
to ~ 10"'^^ G in dark-energy dominated universes [5]. In the absence of dynamo, seeds close to 
10-12 G (or even up to 10 ^ G) are necessary to meet the observations. An additional issue is 
the scale of the initial magnetic field, since the galactic dynamo requires seeds with a minimum 
(comoving) coherence length of ~ 10 Kpc to operate successfully [3j0 

The literature contains various scenarios of cosmic magnetogenesis, which are traditionally 
clasified into those operating after recombination and those advocating a cosmological origin 
for the large-scale magnetic fields of the universe. Primordial magnetism is an attractive idea 
because it makes these large-scale fields (especially those found in distant proto-galactic struc- 
tures) easier to explain. The case was strengthened further when recent observations detected 



^AU scales and strengths have been redshifted to today, unless stated otherwise. 



1 



substantially strong magnetic fields at very high redshifts [6]. Nevertheless, early-universe mag- 
netogenesis is still not problem-free. Producing magnetic seeds that satisfy the above mentioned 
dynamo requirements has, as yet, proved difficult. There are problems both with the scale 
and the strength of the initial field. Causality means that all magnetic seeds generated be- 
tween inflation and (roughly) recombination have coherence lengths well below the minimum 
dynamo requirements. For instance, i?-fields produced during the electroweak phase transition 
span scales of the order of the astronomical unit. A mechanism known as 'inverse cascade' 
can increase the correlation length by transferring magnetic energy to larger scales, but typi- 
cally requires large-amounts of helicity [7j. Inflation has long been considered as an answer to 
the scale problem, since it naturally produces long-wavelength correlations. However, i?-flelds 
that have survived a phase of de Sitter expansion are typically too weak to sustain the galactic 
dynamo, unless classical electromagnetism is modifled [HI [9] , or Friedmann- Robertson- Walker 
(FRW) hosts with non-Euclidean spatial geometry are employed [lOjFI 




Gravitational waves could also provide a geometrical mechanism of magnetic ampliflcation 
on cosmological scales. Gravitons are readily produced in almost all inflationary scenarios and 
interact with magnetic flelds through the shear anisotropy that they induce. The theory behind 
the coupling between the Weyl and the Maxwell fields in cosmology has been analysed in |lHll2j. 
with the latter article arguing for a potentially very strong (resonant-type) amplification of large- 
scale magnetic fields. It should be noted, however that [12j did not explicitly show the amplitude 
of the S-field to diverge, when the latter oscillated in tune with the gravitational waves. Here, we 
take a step further in that direction. We consider two scenarios and provide fully analytical (and 
gauge-invariant) treatments in both cases. The first operates in highly conductive environments 
and the second assumes poor electrical conductivity. In the former case, the gravito-magnetic 
interaction causes the superadiabatic amplification of the -B-field on super-Hubble lengths. The 
standard adiabatic magnetic decay-rate {B oc a^^, where a is the scale factor) slows down. In 
the second scenario, where the conductivity is very low, the Weyl-Maxwell coupling proceeds 
differently. There, the solutions show resonant magnetic amplification on all scales of interest. 
Driven by the gravitational waves, the S-field oscillates with an amplitude that is found to 
diverge when the wavelengths of the two sources coincide. 

The objectives of the present article are (i) to demonstrate that the nature and the out- 
come of the gravito-magnetic interaction depends on the electric properties of the environment, 
since there is considerable confusion on this issue in (some of) the literature, and (ii) to show 
that Weyl curvature distortions can resonantly amplify the S-field. Mathematically speaking 
the latter means that the equations contain resonant solutions. The basic physics behind these 
gravito-magnetic resonances is relatively simple. Gravitational waves are known to trigger shear 
inhomogeneities that affect the magnetic evolution. In poorly conductive environments, where 

^Large-scale magnetic fields in Friedmann universes have long been believed to decay 'adiabatically'. The 
conformal flatness of the FRW models and the conformal invariance of Maxwell's equations are thought to guar- 
antee a Minkowski- like {B oc a~^, with a being the cosmological scale factor) decay-rate for the B-fleld. However, 
although all Friedmann universes are conformally flat, only the spatially flat is globally conformal to Minkowski 
space. For the rest, the conformal mappings are local. As a result, the adiabatic magnetic decay is only guaran- 
teed in FRW models with Euclidean spatial sections. Magnetic fields embedded in Friedmannian universes with 
open spatial sections are 'superadiabatically' amplified on scales close to the curvature radius [10]. The hyperbolic 
geometry slows the decay rate to , leading to astrophysically interesting residual _B-fields without breaking 
away from conventional electromagnetism or abandoning the symmetries of the FRW spacetimes. 




2 



the i3-field naturally oscillates, the aforementioned shear perturbations act as an external driving 
force. Forced vibrations, however, are well known to provide the physical stage for resonances to 
occur. It is not surprising then that there exist analytical solutions showing the magnetic ampli- 
tude to diverge when the Weyl and the Maxwell fields oscillate in tune. Realistically speaking, 
of course, what such diverging solutions indicate is an increased chance of a disproportionately 
strong amplification. This is also what our analysis suggests: that it is theoretically possible to 
substantially amplify cosmological S-fields by means of very weak Weyl-curvature distortions 
alone. 

2 The Weyl-Maxwell coupling in cosmology 

In a general spacetime, with metric gab of signature (— , +, +, +), we define a timelike vector field 
u"" = dx^/dr that is tangent to the worldlines of the fundamental observers (with r representing 
their associated proper time). Together with the projector, hat = 9ab + UaUb, the Ua-field 
introduces a unique 1+3 threading of the spacetime into time and space, while it decomposes 
every variable, every differential operator and every vector/tensor equation into their timelike 
and spacelike components (see [13] for details). 

2.1 Background dynamics 

Consider a spatially flat FRW spacetime filled with a single barotropic fluid and permeated 
by a weak magnetic field (Ba). The latter has energy density well below that of the matter 
(i.e. = BaB"^ <C p) to ensure that the symmetries and the evolution of the of the host universe 
remain unaffected. We may therefore ignore terms of order B^ and the magnetic contribution 
to the background dynamics. The zero-order model is then determined by the standard set 

H^ = Ip, H = -H^-Ip{1 + 3w) and p= -3{l + w)Hp, (1) 
6 

where w = p/p is the barotropic index of the fluid {p is its isotropic pressure), H = a/a is the 
Hubble parameter and overdots indicate (proper) time derivatives!! At the same time, 

4 = -2HBa , (2) 

monitors the zero-order evolution of the (test) B-Held |14j . 

Next, we will perturb the above background by allowing, among others, for gravitational- 
wave distortions. To first order, these are monitored by shear perturbations (see ^ 12.21 below), 
which couple to the magnetic field through the induction equation p[l]. Our aim is to study 
the effects of gravitational waves on the evolution of large-scale magnetic fields during the early 
(pre-recombination) stages of the expansion. The covariant expressions governing the gravito- 
magnetic interaction have been derived in [121, where the reader is referred to for the details. 
Here, we will simply present the relations and use them to study the Weyl-Maxwell coupling in 
cosmology. Further discussion on the formalism and its applications can be found in |14j . 

^Throughout this paper we use geometrised units with SttG = 1 — c. 
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Before closing this section, it is worth noting that we follow the approach of [IH I12j . which 
also allows for a weak background (test) magnetic field. Alternatively, one could treat both the 
-B-field and the gravitational waves as first-order perturbations and study their interaction at 
second order [15 j . The latter scheme maintains the purity of the FRW background, but only 
allows shear-magnetic terms into the equations. All other second-order couplings are excluded. 
As a result, when it comes to the evolution of the i?-field, the two approaches are equivalent. 
The second-order magnetic wave equation of [15] is the same with that of [12] (compare Eq. 
(45) in [15] - see also comment immediately after expression (47) there - with Eq. (5) in |12) . 
or with relations ([5]), pO|) here). 

2.2 The gravito- magnetic interaction 

Information on propagating gravitational waves is encoded in the non-local component of the 
Riemann curvature tensor, which is commonly referred to as the Weyl field. The latter splits into 
an electric and a magnetic part, represented by the symmetric and trace- free tensors Eat and 
Hab respectively. The transverse components of these quantities describe gravitational waves, 
with their transversality guaranteeing that only the pure-tensor modes are accounted for and 
all scalar and vector perturbations have been switched off. In that case the magnetic Weyl 
component is directly related to the shear by means of the linear constraint |13] 

Hab = cuvlaab , (3) 

with cuvlSab = £cd{a^'^Si,^'^ for any orthogonally projected, symmetric and trace-free tensor Sab- 
All these mean that, to first order, we can monitor gravitational waves using the wave equation 
of the shear. On a spatially flat FRW background, the latter linearises to \\.2\ [T3] 

dab - DV^b + mdab + ^ (1 - ^w)H^aab = , (4) 

where = D^D'^ is the orthogonally projected Laplacian operator]^ 

The electromagnetic field obeys Maxwell's equations which combine to provide a set of wave 
equations, one for each component (see expressions (39) and (40) in [16j). After all scalar and 
vector perturbations have been switched off, and given our weakly magnetised background, the 
wave formula of the i?-field reduces to 

Ba - B^Ba + 5HBa + 3(1 - w)H^Ba = 2 i&ab + 2Haab) & + curlj-^ . (5) 

The above monitors the gravito-magnetic interaction on a spatially fiat FRW background at the 
lowest perturbative order [12] . Note that Ba is the background S-field (see § 12. ip and Ja the 
electric 3-current. This depends on the electric properties of the medium, which are encoded in 
Ohm's law. In the frame of the fiuid (i.e. that of the fundamental observers), the latter takes 
the covariant form 

Ja = ^Ea , (6) 

^In [15] the shear was also used to describe gravitational waves at second order. Beyond the linear level, how- 
ever, expression ((3]) no longer holds and the shear alone cannot (generally) monitor gravitational- wave distortions. 
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where is the electric conductivity of the matter [T7]. In poorly conductive environments, <; — )• 
and there are no currents. Then, we may ignore the last term in the right-hand side of ([5]). At 
the opposite end, where the conductivity is very high and ? — )• oo, the electric field vanishes. In 
that case, J a = cuilBa (e.g. see [H]) and expression ^ reduces to 



Not surprisingly, in the absence of its electric counterpart, the 5-field no longer obeys a wave-like 
formula. Instead, the magnetic evolution is determined by the familiar induction equation |14] . 

3 Gravitationally driven magnetic field 

The significant differences between expressions ([5]) and ([7]) suggest that the electric properties of 
the cosmic medium, at the time and in the region the gravito-magnetic interaction occurs, are 
crucial. Here, we will mainly focus on the two limiting cases of very low and very high electrical 
conductivity. The intermediate case of finite resistivity will also be discussed in § 14. 3i 

3.1 Harmonic splitting 

Gravitational waves introduce shear inhomogeneities, which couple to the background magnetic 
field through Maxwell's formulae. The resulting gravito-magnetic stresses drive the evolution of 
the perturbed (the induced) i?-field according to Eqs. ([5]) or ([7]), depending on the conductivity 
of the environment. Our next step is to harmonically decompose these expressions. This will 
also provide a relation between the wavelengths of the agents involved. Splitting the zero-order 
field as Ba = i?(„)Qa"\ helps to assign a scale and a finite wavenumber (n) to it. Similarly, 

(k) (£) 

we may set aab = '^k(^{k)Qab ^^^'^ ~ '^iB(£-jQa , where k, I are the associated wavenumbers 
and Q^a^ = Q^^J Q^^^ by construction [TTl [12] . The three corresponding wavevectors satisfy the 
relation = = (n^ + ka){n'^ + k°'). This immediately translates into 



with < ip < IT representing the angle between the wavevector of the background magnetic field 
and that of the oncoming gravitational wave. 

3.2 Evolution equations 

The harmonic decomposition of the previous section reduces the wave formulae and ([5]) 
from partial to ordinary differential equations. Together with the decomposed counterpart of 
expression ([7]), the new set consists of 



Ba + 2HB, 



(TabB 



(7) 




(8) 




(9) 
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and 



3{l-w)H'+ - 



Br 



2 + 2Ha 



{k)) ^(n) 



(10) 



+ '^HB^^g) = 5(„)(T(fc) . (11) 

The first two monitor the gravito-magnetic interaction in poorly conductive cosmological envi- 
ronments, while the set Q, (jlip does the same when the electric conductivity is high and the 
ideal magnetohydrodynamic (MHD) limit applies. 

We proceed further by introducing the dimensionless, expansion normalised variables ^{k) = 
a(^i;-^/3H and = B(^(^/2,H. Then, using conformal (r/, with i] = 1/a) instead of proper time, 
the system (|9]l- (fTT]l recasts into 



^k) + (1 



3w) 



3[1 + (2 - 2,w)w] 



(k) 



0. 



(12) 



3(1 — 3w;)tt; / a' 
a 



B 



2alB, 



in) 




and 



B' 



1 



W + 2 



3^)(^ 



^\k) + 



1 -3'u; 



B, 



(^) 



^(fc) , 



(13) 



(14) 



where primes indicate conformal-time derivatives and the zero suffix marks the beginning of the 
gravito-magnetic interaction. Also, in deriving Eqs. (jl3l) and (I14p we have used the evolution 
law, B(^n) ^ of the background magnetic field (see Eq. ([2]) in ^ 12. ip . As before, the set (fT2]) . 



(fT3l) applies to poorly conductive environments, while Eqs. (fT2]) and (fT4l) operate at the ideal- 
MHD limit. The advantage of the new notation is that, when evaluated in the radiation era, 
expressions (I12p - (114p take a very simple and compact form that is also straightforward to solve 
analytically. In addition, during the radiation era, the expansion normalised i3-field satisfies the 
criteria for gauge invariance (see § H] below) . 

Mathematically speaking, expression (I13p is a wave formula and (I14p is a first-order differen- 
tial equation. Physically, this difference reflects the absence of electric fields in highly conductive 
environments. What is most important, however, is that in (I13p - as well as in (llOp - the gravi- 
tational waves drive the magnetic oscillation. Forced oscillations, however, are known to provide 
the natural physical stage for resonances to occur. In other words, expressions (llOp and ()13p 
open the theoretical possibility of gravito-magnetic resonances on cosmological scales. 



4 Gravitationally amplified magnetic field 

The galactic dynamo requires seed fields coherent on a minimum comoving scale of roughly 
10 Kpc. In typical cosmological models such scales have been outside the Hubble horizon for a 
long time. Therefore, in order to achieve astrophysically interesting results, we should consider 
lengths close and beyond the Hubble radius. 
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4.1 Resonant amplification 

Inflation has long been known to provide a mechanism that produces electromagnetic fields and 
gravitational waves with super-horizon correlations. Also, throughout the inflationary phase, 
the universe is a very poor electrical conductor. Although the conductivity is expected to grow 
quickly with the onset of the radiation era, causality must confine its efi^ects (i.e. the currents 
that will eliminate the electric fields) well inside the horizon. Near and beyond the Hubble 
radius, the gravito- magnetic interaction will still proceed in line with Eqs. (|12p and (jl3p . 

In the radiation era the expansion normalised variables introduced in § 13.21 have the dou- 
ble advantage of simplifying the equations and also of freeing the results from gauge-related 
ambiguities. The gauge issue concerns the magnetic field, since the shear vanishes at the zero- 
order level and all the shear-related variables are therefore gauge independent. The presence 
of a nonzero background magnetic field, on the other hand, means that the magnetic-related 
variables are generally gauge dependent. During the radiation era, however, the scalar i?(£)/3i? 
remains constant to zero order and therefore satisfies the Stewart &: Walker criteria for gauge 
invariance [18]. As a result, when radiation dominates the energy density of the universe, the 
expansion normalised variable B(^£-j is a gauge-invariant linear perturbation. 

With these in mind, we apply the system ()12p . ()13p to a radiation dominated almost-FRW 
universe. Setting w = 1/3, a = aQHot] and a' /a = 1/r], Eq. (fT2|) reduces to 



(15) 



and accepts the oscillatory solution 



= Ci 



sm{kr]) + 



cos{kr]) 
krj 



+ 0, 



cos(A;r7) 



krj 



(16) 



Keeping only the dominant S-mode, we may rewrite the above a^ 



-(fc) 



sin(fcr/o) + ^ Sq^'"'^ cos{kr]o) 



+ 



cos(A;r/o) - ^ Sq^''^ sin(A;77o) 



sin{kri) 
cos{kr]) . 



(17) 



with the zero suffix indicating the onset of the gravito-magnetic interaction - which for our 
purposes coincides with the beginning of the radiation era. 

At the same time and on the same scales, the gravitationally driven magnetic field obeys the 
wave formula 

r/ 



(18) 



with B(^n) — B{n)/'^H and "^(k) given by solution ()17p . Differentiating the latter with respect to 



>(n)/ 



^Including the subdominant E-mode in (|17[) simply adds decaying oscillatory terms to the right-hand side of 
solutions (|19l) and (|20|) . without changing neither the nature of these expressions nor the overall results. 
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conformal time and substituting the result into the right-hand side of (jlSp . we arrive at 



.55 s, 

rj 



(n) 



T.f'' cos{kr]o) - i Sq^''^ sm{kr]o) 



Sq''^ s\a{kr]Q) + ^ Sq'^"-' cos(/cr/o) 



1 



v(fc) 



sin(/c7/) 
cos{krj) . 



(19) 



The above traces the evolution of the expansion normalised, gravitationally driven magnetic field 
during the radiation epoch of an almost-FRW universe near and outside the Hubble horizon. 
The solution of Eq. leadE 



71 



■ Si[(.^ T k)rj\ sin [{Irj =F krjo) =F 



n 



i 



■ Ci[{£ =F A:)/?] cos [{iri =F /c??o) =F 



(20) 



where tan( 



/'*^^o 



||see Appendix for the details). Also, Si(3;), Ci(x) represent the sine 
and the cosine integral functions respectively [19]. Therefore, driven by the gravitational waves, 
the expansion normalised, gravitationally driven magnetic field oscillates with an amplitude that 
depends on the initial conditions. As typical for forced oscillations, the amplitude also depends 
on the wavenumbers of the agents involved. Given that 



lim Si(x) 







and 



lim Ci(x) 



-oo . 



(21) 



we deduce that when the two waves are in resonance (namely as k ^ i) the amplitude of the 
magnetic field increases arbitrarily. Note that, at least theoretically speaking, resonant effects do 
not depend on the strength of the driving source. In this respect, the magnetic amplification seen 
in solution (|20p requires only the mere presence of gravitational waves and is rather insensitive 
to the amount of energy stored in them. 

We may also express solution (1201) in terms of the actual (instead of the expansion normalised) 
variables and with respect to proper (instead of conformal) time. To do that recall that = 
S(£)/3i?, E(fc) = (T(fc)/3i/ and B(^n) = .S(n)/3^ by definition, while S'^^-j = aa(^i^^/3H after keeping 
only the dominant shear mode. Also, H = l/2t and r/ = \/i/aoHQ^/to during the radiation era. 
Then, after some lengthy but straightforward algebra, we arrive at 



B. 



±- 



B, 



in) 



B, 



in) 



Si 

\aoHo 




to. 



sm 



£Ho 



tj \aoHoy toj 



cos 



t k 



(22) 



where now tan(/> = aa^' /ka^\ The above shows that the gravitationally driven (actual) 



magnetic field decays adiabatically (i.e. oc t 



oc a 



like its background counterpart. 



^Mathematically speaking, Eqs. p9|) and (|20|l hold on all scales, although in this section the 'operational' 
domain of these expressions lies near and outside the Hubble radius. This does not always need to be the case 
however. As we will see in § 14.31 Eq. (|19|) and solution (|20|) can also extend to subhorizon lengths. 
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However, as in solution (j20p . the presence of the cosine integral function Ci[(^ — k)\/t/aQHQ^/tQ\ 
in the right-hand side of (j22p guarantees that the amplitude of the S-field diverges when the two 
sources are in resonance (i.e. for k ^ £). It goes without saying that solution (|22|) can be obtained 
directly from the system ([9]), (jlOp . In particular, solving Eq. keeping the dominant shear 
mode and then substituting into pop leads to expression (|22p . To the best of our knowledge, 
this is the first time the amplitude of cosmological (or astrophysical) magnetic fields has been 
found to diverge, as a result of their interaction with gravitational-wave distortions. 

Solutions (|20p and (j22p describe one type of resonance between the Weyl and the Maxwell 
fields, but there are other possibilities as well. In each individual case the characteristics of the 
resonance are decided by the shape of Eq. (jlSp and by the form of the driving term in the right- 
hand side of that expression. Suppose that the source term has time-independent amplitude. 
This can happen, for example, when the background magnetic field decays as r]~^ (like many of 
the -B-fields proposed in [8]- [10]), while the expansion normalised shear is still given by ()17p . In 
that case Eq. (fT8]l accepts the solution 



which also ensures arbitrary growth for the amplitude of the gravitationally driven i3-field as 
k ^ i. It is therefore plausible to argue that a greater variety of physical environments could 
support analogous effects. It is also very likely, however, that only a small fraction of these cases 
will be analytically tractable. Overall, our analysis opens the theoretical possibility for gravity 
to act, through its long-range component, as a very efficient magnetic dynamo on cosmological 
scales and long before the onset of structure formation. 

In retrospect, the results reported in this section should not come as a surprise. After all, 
we were dealing with forced oscillations, where the Weyl part of the gravitational field was the 
driving agent. Forced oscillations have long been known to provide the physical stage where 
resonances naturally occur. In this sense, our study simply shows that resonances can also take 
place between cosmological gravitational waves and large-scale magnetic fields. 

Finally, before closing this section, we should underline that solution (j22p and the resulting 
gravito-magnetic resonance could have been obtained in [15] as well. For example, solve the 
linearised version of Eq. (46) - the numbering is that of [15] - during the radiation era and on 
all scales. Then, substitute the wave solution of (46) into Eq. (45) - the numbering is still that 
of [15] - and solve again. The result would have led to expression (j22p above. 

4.2 Superadiabatic amplification 

Let us now turn to the ideal-MHD limit. The main objective is to illustrate the qualitative change 
in the nature of the gravito-magnetic interaction caused by the different electric properties of the 
cosmic medium. To some extent this is to be expected, given the significant differences between 
Eqs. dSj) and ^ - or between expressions (fT3]) and ([H]). To keep things simple, we assume that 
the magnetic field is frozen-in with the matter on all scales. This is possible within, say, the 
Eddington-Lemaitre or the 'emergent-universe' scenarios [20], where the cosmos evolves from 
a past-eternal static state. Then, in the frame of the fluid, the gravito-magnetic interaction is 




(23) 
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monitored by Eqs. (jl2p and (jl4p . This means that during the radiation era the gravitationally 
driven magnetic field evolves according to 

B[,^ = Uj-%,y (24) 

with = ^(-„)/3i? and Ej-;,) given by ^TU\i . As before, the zero suffix indicates the beginning 
of the interaction. The form of expression (j24p leads to the first qualitative conclusion. In 
particular, the fact that the -B-field no longer obeys a wave equation means that gravito-magnetic 
resonances are not possible in highly conductive environments (see also § VIB in [12j). 

To demonstrate the effect of the Weyl field at the ideal MHD limit, it suffices to consider 
super-Hubble lengths and replace solution (jl6p with its long-wavelength approximation, namely 




Substituting this into Eq. (j24[) . integrating and setting the integration constant to zero for 
simplicity, leads to 

s„ = isr(4''+')»4''')(^)'-fii"'(24''-*4''')(*) . 

on all super-Hubble lengths. Finally, recalling that = i?(^)/3i7, that = i3(„)/3i/ and 
that H = l/cgi^o??^ throughout the radiation era - which also implies that ao-ffo% = 1 at the 
time, we arrive at 

i^W = \ {^t' + m^r) - (24'=^ - Vo^r) (^) ' • (27) 

Following ()27p . prior to equipartition, the gravitationally driven magnetic field no longer 
decays adiabatically. Instead, the dominant magnetic mode remains constant, which impies 
a superadiabatic type of amplification for the i?-field on all super-Hubble scales due to Weyl 
curvature effects. This means that magnetic fields that cross inside the Hubble radius around 
the time of matter-radiation equilibrium have magnitudes equal to 

B^§ = \{^t^+n,^r)Bt'- (28) 

Thus, the strength of the residual i?-field depends on that of its background counterpart at the 
onset of the interaction and on the relative shear anisotropy at the same time. In particular, 
the final outcome is decided by the values of t!"^^ and of which the first (at least) is 

expected to be very small. Therefore, unless ryo^Q » Sq - or the shear anisotropy is larger than 
anticipated - this effect is unlikely to produce magnetic fields of astrophysical importance. 

4.3 Finite conductivity 

So far, we have focussed on the two limiting cases of very low and very high electrical conduc- 
tivity. For completeness, let us take a brief qualitatively look at the intermediate scenario of 
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finite resistivity. Assuming that the conductivity coefficient has a nearly homogeneous spatial 
distribution (i.e. setting Dq? ~ 0), Eqs. ([5]) and (j6]) combine to 



Ba-B'Ba + 5 ( 1 + ^ ? j HBa + 3 (l-w + ^^] H'Ba = 2 



&ab + 2[l + ^^] Ha ah 



B\ (29) 



where <; = q/H is the expansion normalised conductivity scalar (see § VIA in |12] for more details 
on the derivation). As long as ? ^ 1, the above reduces to Eq. (jlOp and the Weyl-Maxwell 
coupling proceeds as in poorly conductive mediums irrespective of the scale. Therefore, gravito- 
magnetic resonances, like those discussed in § 14.11 before, can also occur on subhorizon lengths, 
as long as <; ^ 1 there 

When the ratio q/H \s large, we can no longer ignore the role and the effects of the electric 
conductivity (on sub-Hubble scales). Then, expression (|29p assumes the form 

Ba - B^Ba + ^HBa + 2^H^Ba = 2 {dab + \ ^Haab^ & , (30) 

and the gravito-magnetic interaction depends largely on the electric properties of the medium. 
Situations like this require a model for the conductivity of the cosmic plasma. Note that the 
above relation applies to cosmological environments of significant but still finite conductivity 
and one should not simply extrapolate it to the ideal-MHD limit (where — )• oo). There, the 
electric fields vanish (see Eq. (^) and one should instead use Eqs. PT]) or ([HD. 



5 Discussion 

The widespread presence of magnetic fields in the universe is a hard observational fact. Prom 
the nearby planets and stars, all the way to distant galaxies, galaxy clusters and remote high- 
redshift proto-galactic clouds, the magnetic presence has been repeatedly verified. The origin of 
the large-scale S-fields, however, is still illusive. Typically, magnetic fields with coherence scales 
large enough to sustain the galactic dynamo are too weak, while those of adequate strength 
span very small lengths. Inflation can naturally produce large-scale correlations and it has 
long been seen as the strongest candidate for early-time magnetogenesis. Nevertheless, -B-fields 
that survive the exponential phase are typically too weak to seed the dynamo, unless classical 
electromagnetism is modified or FRW hosts with non-Euclidean spatial geometry are involved. 

Geometry could provide a solution to the magnetic-strength problem through Weyl-curvature 
distortions as well. Gravitational waves, which are readily produced during the inflationary era, 
interact with magnetic fields via the shear anisotropy that they introduce. The outcome of 
the coupling between the Weyl and the Maxwell fields depends, among others, on the electric 
properties of the medium that fills the universe. In highly conductive environments, where the 
currents keep the S-field frozen-in with the fluid, the gravito-magnetic interaction slows down 
the standard adiabatic magnetic decay. Looking at the radiation era, we found a superadiabatic- 
type of amplification. The strength of the field remained constant, instead of dropping in accord 
with the usual a~^-law. Nevertheless, the overall magnetic growth is likely to be small, unless 
the amount of the initial shear anisotropy is higher than expected. 
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The Weyl-Maxwell coupling looks much more promising in poorly conductive environments, 
where the gravitational waves drive the magnetic oscillations. Forced oscillations are well known 
to provide the natural environment for resonances to occur. It is not surprising then that our 
analysis has found resonances between gravitational waves and magnetic fields. During the 
radiation era, in particular, the amplitude of the gravitationally driven i?-field was found to 
diverge when the wavenumbers of the two sources coincided. Physically speaking, what resonant 
solutions mean is that very weak vibrations can have disproportionately strong effects. This 
is also what our analysis indicates: that very weak gravitational waves can strongly amplify 
cosmological magnetic fields. 

Typically, the mechanism operates near and beyond the Hubble horizon. Smaller scales can 
also provide a favourable environment, as long as? = ?/-ff<Cl there (see Eq. (I29p in § 14. 3p . 
Given that inflation produces gravitational waves and magnetic fields essentially on all lengths, 
the chances of resonances to occur (i.e. the probability of their wavelengths coinciding) should 
be relatively high. On causality grounds, the most likely time for the amplification to occur is at 
horizon crossing. Well inside the Hubble radius, the currents should eliminate the electric fields 
and freeze the magnetic component in, thus changing the nature of the interaction. Naively, one 
should also expect a similar amount of strengthening for every magnetic mode passing through 
the horizon. This favours a scale-invariant spectrum for essentially all the large-scale i?-fields 
in the universe. The latter appears to be in qualitative agreement with the observations, which 
show magnetic fields of similar (//G-order) strength in nearby galaxies as well as in remote 
proto-galactic structures. 

Finally, we should point out that cosmology is probably not the only venue where gravito- 
magnetic resonances can occur. In principle, analogous effects could also take place in astrophys- 
ical environments. After all, the nature and the basic physics of the interaction remain the same. 
With this in mind, it might be worth investigating whether (and under what conditions) the 
coupling between astrophysical gravitational waves and magnetic/electromagnetic fields could 
lead to the resonant amplification of the latter. 

The author thanks Axel Brandenburg, Ruth Durrer and Elisa Fenu for helpful discussions. 

Appendix 

Given the key role of solution ()20p , we provide some additional technical information and inter- 
mediate steps for the interested reader who wishes to recover the above named expression. To 
begin with, it is relatively straightforward to verify that Eq. (jl9p can be solved analytically on 
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all scales. The solution describes forced oscillation and takes the form 



t3[e) = Ci sm{ir]) + C2 cos{ir]) 
+ 



^- — — [cos{kr]o) sm{ir]) — sin(A;?7o) cos(£r/)] Si [{i — k)r]] 



^ ^ ^ — [sin(A;r/o) sm{£r]) + cos(A;?7o) cos(£ry)] Si [(£ — k)r]] 

H ^ — — [sin(A;77o) sin(£r/) + cos^krjo) cos(£r/)] Ci [(^ — k)ri] 

H ^ £ ^ — [cos(A;r/o) sin(£r/) — sin(A;?7o) cos(£?7)] Ci [{£ — k)r]] 

^ — — [cos(A:t/o) sm{£r]) + sin(A:T/o) cos(£rj)] Si [(^ + /i;)ry] 

— ^ — [sin(/i;r/o) sm{ir]) — cos(A;?7o) cos(£ry)] Si [{£ + k)r]] 
^ — — [sin(/c7/o) sin(^r/) — cos{kr]Q) cos{£r])] Ci [(^ + k)r]] 



+ 



+ 



H — [cos(A;r/o) sin(£r/) + sin(A;?7o) cos(£r?)] Ci [{£ + k)r]] 



(31) 



where Ci^2 are the integration constants, while Si(x) and Ci(x) are the sine and the cosine integral 
functions respectively. Also, primes indicate conformal time derivatives and the zero suffix 
corresponds to the onset of the gravito-magnetic interaction, which for our purposes coincides 
with the beginning of the radiation era. Since the gravitationally induced magnetic field vanishes 
in the absence of shear perturbations, we may set the aforementioned integration constants to 
zero. Then, using standard trigonometry, solution ([3T]) reduces to the compact expression 



3B, 



5(n) 



kY^-Q^ sm.{lri =F fcr/o) =F '^'q'^' cos{£r] =F k-qo) Si [{£ =F k)ri] 



±- 



3B, 



(n) 



kJ:l^^ cos{£ri =f kr]o) ± Sq^"^ sin(^7? =f krjo) Ci [{£ ^ k)v] 



v(fc) 



(32) 



f (k) (k) 

Finally, setting tan0 = Sq /kT,Q and employing some relatively lengthy but straightforward 



algebra, we arrive at solution (pOj) 



B 



(^) 



£ 



Si[{£ =F k)7]] sin [{£r] =F kijo) =F 



Ci[{£ =F k)ri] cos [{£rj =F kr]o) =F 



(33) 



Note that, for simplicity, we have suppressed the wavenumbers (k) of the Sq modes inside the 
square roots. 
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